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We derive the multi-fractal scaling of probability distributions of multi-particle configurations for 
the binary reaction-diffusion system A + A — ► in d < 2 and for the ternary system 3A — > in 
d — 1. For the binary reaction we find that the probability Pt(N, AV) of finding TV particles in a 

fixed volume element AV at time t decays in the limit of large time as (^) (hit) 5 for d = 2 

and t ' t s y ' for d < 2. Here e = 2 — d. For the ternary reaction in one dimension we 

i i jv/o N{N — 1)(JV — 2) 

find that Pt(N, AV) ~ (i^)™ /2 (ln i) 3 . The principal tool of our study is the dynamical 

renormalization group. We compare predictions of e-expansions for Pt (TV, AV ) for binary reaction 
in one dimension against exact known results. We conclude that the e-corrections of order two and 
higher are absent in the above answer for Pt(N, AV) for N = 1, 2, 3, 4. Furthermore we conjecture 
the absence of e 2 -corrections for all values of N. 



I. INTRODUCTION 

Modeling chemical reactions is an important practical 
and theoretical problem. Systems of reacting particles 
are typical of complex irreversible non-equilibrium sys- 
tems. The popular description of these systems in terms 
of simple rate equations currently adopted in chemical 
kinetics [lj does not appear to have a firm theoretic foun- 
dation and often produces wrong results. 

Interacting particle systems provide a good model for 
simple chemical reactions. Well-known examples include 
systems of diffusing-coalescing and diffusing-annihilating 
particles describing reactions A + A — ► A and A + A — > 
respectively. Extensive studies of these particle systems 
in low dimensions have shown that the rate equations 
yield incorrect results for the computation of average 
concentrations of reactants, see Q for review. Rate 
equations fail in low dimensions due to the presence of 
large fluctuation effects, which violate mean-field theory 
(MFT) assumptions underlying their derivation. The 
study of diffusive annihilation (or coalescence) ^{J is a 
good starting point for analyzing large fluctuation effects 
in more complicated non-equilibrium statistical systems. 
For example, the system A + A — > A was used to analyze 
the aggregation of massive diffusing particles [f| . 

The aim of the current paper is to study the effects of 
large fluctuations on correlation functions of an arbitrary 
order for the reactions A + A— >• in <i < 2 and A + A + 
A -> in d = 1. 

Large fluctuation effects are accounted for in binary 



reaction-diffusion models using Emptylnterval methods 
(EIM) and its generalizations |g, 0, H Q ■ This approach 
is restricted to d = 1 and does not extend to higher 
dimensions. There are rigorous results on the average 
density of particles in d = 1,2 . The Smoluchowski 
approximation gives correct answers for average concen- 
trations but cannot be used for higher order correla- 
tion functions. In the early 90's, the work of Cardy and 
Lee used field theoretic methods, in particular the 

renormalization group (RG) to obtain an answer for the 
average density as well as its amplitude for d < 2. The 
study of A + A — ► 0(A) has also introduced the concept 
of stochastic rate equations with imaginary multiplicative 
noise 0. 

In this paper we derive the multi-scaling of correlation 
functions in the systems A + A — > and 3A — > using 
the RG method. The main object of our study is the 
large time temporal scaling of the probability of finding 
N particles in a small volume AV, denoted Pt(N, AV). 

For the A + A — > reaction-diffusion system we are 
interested in d < 2. We do not consider higher dimen- 
sions as the answers there are given by MFT. Most stud- 
ies have concentrated on computing the average density 
of particles (N = H^. To the best of our knowl- 
edge, the computation of multi-particle probabilities are 
only considered in 0, Q, 13 > with the analysis restricted 
to one dimension. Dynamical RG method allows us to 
obtain answers for the large time limit in the form of 
an e-expansion (e = 2 — d) for d < 2 and logarithmic 
corrections to the MF scaling for d = 2: 
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Equation £Q represents the multi-scaling or the deviation 
oiP t (N,AV) from P t (l, AV) N . As 
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equation Q reflects the anti-correlation between parti- 
cles in the large time limit. 

For the ternary reaction-diffusion system 3^4 — > we 
restrict our attention to d = 1. MFT provides the an- 
swers in higher dimensions ^3 • The average density has 
been studied in 0, j and the two-point function in 
fl2j | . Higher order correlations have not been considered. 
The method of empty intervals does not extend to the 
ternary reaction. The RG method produces asymptoti- 
cally exact results for d = 1: 
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The paper is organized as follows: an introduction to 
the lattice model of A + A — > is given in Section II. 
A, its field-theoretic re-formulation is given in Section 
II. B, followed by a summary of the mean field results in 
Section II. C. Section III contains the RG analysis of the 
binary annihilation system for d < 2, starting with the 
description of renormalization procedure in Section III. A. 
Section III.B contains the derivation of the corresponding 
Callan-Symanzik equation for the theory, its solution and 
the derivation of large time asymptotics of multi-particle 
probabilities. We prove the exactness of our result for 
Pt(N — 2, AV) using the first Hopf equation of the the- 
ory. In Section III.C we compare our results in d = 1 
against results from 6] and establish their equivalence 
for N = 1, 2, 3, 4. Further we conjecture the exactness of 
our one-loop answer in d = 1 for general values of N. In 
Section IV we extend our analysis to the ternary reaction 
in d = 1 using the same methodology. 



II. FIELD-THEORETIC FORMULATION AND 
MEAN-FIELD LIMIT OF A + A -> MODEL 

A. The Model 

Consider a set of point particles performing random 
walks, characterized by diffusion coefficient D, on the lat- 
tice Z d . Any two particles positioned at the same site can 
annihilate each other according to an exponential process 



with rate A. For the simplicity of our analysis we assume 
finite reaction rates. However the large time asymptotics 
of our model belongs to the universality class of instanta- 
neous annihilation-diffusion model (see end of Section III 
B). It is assumed that the initial distribution of particle 
number at each site is independent Poisson with mean 
N a . 

Let the random variable N t (x) represent the occupa- 
tion number for site x at time t. The configuration vector 
N_ = {N(x)} xeZ d specifies the state of the system at time 
t by encoding the occupation number at all sites. N_ is 
also termed a microstate. Let Vt(N) be the probability 
of finding the system in microstate N_ at time t. Corre- 
lation functions of Nt(x) can be obtained by averaging 
functions of N(x) with respect to Vt (N). For example, 
the average density is given by: 



(4) 



Due to translational invariance this will be independent 
of x. Our main object of interest is the large time limit 
asymptotics of Vt(N(x) = N), the probability of finding 
N particles at site x. It turns out that in the low density 
limit this probability is proportional to the N th factorial 
moment of N t (x), which we denote by M^(x,t). This 
can be verified as follows: 



M N (x,t) 
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In the above derivation, the first five lines are exact re- 
lations. The last line is due to the fact that in the large 
time limit we expect the particle density to be low and 
particles are anti-correlated 0, 0, B 113 • Hence the con- 
figurations with the smallest possible value of N t = N 
will give the dominant contribution. 

The coarse-grained counterpart of V t [N(x) = N] is 
P t (N, AV), the probability of finding N particles in the 
volume element Ay. Let AN t (x) be the number of par- 
ticles in a volume AV (centred at x) at time t: 
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where nt(y) stands for the density of particles at time 
t at a point y. It should be mentioned that upon the 
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averaging AN t is independent of x due to translational 
invariance. In the limit of large time and fixed AV , fac- 
torial moments of AN are related to Pt(N, AV) via a 
relation analogous to JSJ: 



Let us introduce the quantity A(j>: 



P t (N, AV) 



N\ 
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As we will demonstrate in the next section, the facto- 
rial moments of AN t (x) admit a simple representation in 
terms of polynomial moments of Doi's fields. 



B. Path-Integral representation 

In the last section it was shown that to obtain corre- 
lation functions of interest, we need to know moments 
of N(x) with respect to Vt(N). It is possible to evalu- 
ate these moments in a path-integral setting. The for- 
malism is due to Doi-Zcldovich-Ovchinnikov and we re- 
fer the reader to Q for details. We simply present a 
schematic derivation. The time evolution of the prob- 
ability measure Vt on the space of microstates is given 
by the master equation. The master equation is a linear 
first order autonomous differential equation with respect 
to time, which implies that its evolution operator can be 
expressed as a path-integral. Thus we can find a path- 
integral representation for any correlation function. As 
we are interested in universal properties of the system at 
scales much larger than the lattice spacing, it is conve- 
nient to work with the continuum (coarse-grained) limit 
of the lattice model, which corresponds to some effective 
field theory in imaginary time. The continuum limit in 
the path-integral is taken according to the rules 



\{Ax) d -> A D(Axf 



D 



where Ace denotes the lattice spacing. The reader is asked 
to refer to 12] for more details. The resulting field theory 
is given by the effective action S: 

S = J dtj d d x^{d t ~A)(t) + 2\^cj) 2 + \^ 2 <i) 2 ~n ^8{T) 

where n is the initial average particle density. We will 
be working in the large n Q limit. The diffusion coefficient 
D can be set to 1 as it is not renormalized by fluctuations 
[l2|. The </>-field is a response field. The (/>-field is related 
to the local density field, n t (x) in the sense that there is a 
one-to-one correspondence between correlation functions 
of n t and </>. Let (•) be the averaging with respect to the 
path-integral measure e~ s . Then the average density is 

m 



n t (x) = (4>(x,t)) = V(j)(x,t)V4>(x,t) 4>(x,t)e 
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where the volume AV is centred at x. Averaging i|l(J|) 
results in an x-independent quantity due to the trans- 
lational invariance of the system. Let us consider the 
statistics of AN t (x) defined in equation ©. Moments of 
AN and moments of Acf> are related as follows 



(A0) = AN t 
<A0 2 ) 



ANt(ANt-l) 



(11) 
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The relations JTTJ, O follow from 0. Let us note that 
(|llfl gives us a path-integral representation of the average 
number of particles in a volume AV. More generally |l4| 



<A<^) = M N (t) 



(13) 



where AIn is defined by equation 10. In the previous 
section it was shown that in the large-time limit the fac- 
torial moment is approximately equal to the probability 
of finding N particles in a volume element AV. Hence, 



P t (N,AV) = ±(A4> N ) 



(14) 
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As we mentioned before, we are interested in the limit 
AV — > 0. Physically this limit corresponds to AV <C l d , 
where the correlation length I ~ \Ji 4]. By working in 
the large time limit we are effectively dealing with the 
case AV — > 0. Then we may write: 

(A<f> N ) = / dxx.-.dxN {(p(x+Xi)...(j)(x+x N )) (15) 

J AV 

where 

((j){x+Xi)...<i>(x+x N )) = F(xi . . .x N )((/) N (x)) R ■ 

(l + 0(AV 1/d )) (16) 

F(x 1 . . . xn) is the leading coefficient of Wilson's opera- 
tor product expansion (OPE) The limit AV — > 
corresponds to the ultraviolet (x^s — > 0) asymptotics of 
(<f)(x + xi) . . . 4>{x + xn)), which is given by the renormal- 
ized average of the composite operator <p N . We conclude 
that multi-scaling of Pt (N, AV) is described by anoma- 
lous dimensions of the composite operators <f> . The 
spectrum of anomalous dimensions of cf> N will be com- 
puted in Section III. 

C. Mean Field Limit and Loop expansion 

The Feynman rules for the perturbative computation 
of correlation functions are derived from the action 5, 
equation ((HJ . They are given by Fig ^ The propaga- 
tor G Q is the Green's function of the standard diffusion 
equation. The perturbative expansion of (4> m ) in powers 



of A is given by the sum of all diagrams with m outgoing 
lines; hence diagrams contributing to the mean density, 
(4>), have one outgoing line. The action S must be di- 
mensionless; thus in length (L) units, the dimensions of 
the relevant parameters must be 

[t] = L 2 [4>]=L° m=L- d [\}=L d -\ (17) 
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The critical dimension is d c = 2, where the reaction 
rate is dimensionless. For d < 2 the field theory is 
super-renormalizable where all Feynman integrals con- 
verge. For d = 2 we use dimensional regularization to 
ensure convergence of Feynman integrals. 



FIG. 2: a) The classical density; b) The classical response 
function 



lines ((4> m )) and n loops is proportional to n™-g™. Then, 




FIG. 1: Feynman rules for A + A 




n2(i)(l + E<w£(*) 



(20) 



-field theory 



The bare dimensionless reaction rate is given by 
g (t) = A^ 1_d / 2 , which grows with time in d < 2. A com- 
binatorial argument shows that an n-loop diagram con- 
tributing to the mean density is proportional to 3™ _1 (t) 
|12| : thus in the weak coupling regime the main contri- 
bution to the mean density comes from the sum of tree 
diagrams. This is equivalent to the mean-field approxi- 
mation. This sum (tree diagrams) is also termed the clas- 
sical density, denoted n c i(t). In d < 2 MFT is valid for 
small times given initial density is large. This agrees with 
our intuition: for small times local fluctuations around 
the large mean value of the density are small. At large 
times, g (t) grows and MFT breaks down. To compute 
corrections to the mean-field answers, one must com- 
pute higher loop contributions. In particular this involves 
summing infinite sets of diagrams for a fixed number of 
loops. These can be re-summed in a more compact form 
using the classical density (n c i) and the classical response 
function (G c i). The classical response function consists 
of the sum of all tree diagrams with one outgoing and 
one incoming line. The diagrammatic form of the inte- 
gral equations satisfied by these two quantities are given 
in Fig |21 The solution to these equations are: 



n c i{i) 

G c i(x 2 ,t 2 ;xi,ti) 



1 

n c i{t2) 



n c i(h 



(18) 



G (x2,te;xi,ti)(19) 



Using n c i and G c i in combination with vertices of Fig 
1, we arrive at Feynman rules generating finitely many 
diagrams for a given number of loops. 

We use n c i to write a mean-field expansion for (<p m ) 
which is valid for small times in d < 2. Using dimensional 
analysis it can be verified that a diagram with m outgoing 



In this form we see that for small g {t) the loop correc- 
tions are small. Then we may formulate the following 
mean-field answer for the probability of finding N parti- 
cles in volume AV: 

Pt(N,AV) = ^y(A^) 
NV Vl 



N.-N 



(21) 



Comparing (|21(1 with results from 0, Q m d — 1 we 
find the linear temporal scaling of MFT to be incorrect in 
the large time limit. We will compute the correct scaling 
in the subsequent Section. 



III. RENORMALIZATION GROUP ANALYSIS 

OF P t (N,AV) 

The dynamical renormalization group method allows 
one to extract large time asymptotics of correlation func- 
tions of the theory JHJ. The first step is to eliminate all 
e — > singularities of Feynman integrals at some refer- 
ence time t a . The process of removing these divergences 
is called renormalization. This is done by introducing a 
renormalized reaction rate, renormalized fields, etc. The 
number of renormalization constants needed to eliminate 
all divergences is finite for d < d c due to renormalizability 
of (JSJl . Individual terms in the renormalized perturbative 
series for any correlation function C(t) depend on the ref- 
erence time t . The lack of dependence of the unrenor- 
malized version of C{t) on the unphysical parameter t a 
leads to renormalization group (Callan-Symanzik) equa- 
tion for the correlation function. This is solved subject 
to the initial condition at t Q given by the perturbative 
expansion for C(t) at t = t a . This procedure is equiva- 
lent to re-summing all leading e-s ing ularities in d < d c at 
all orders of the loop expansion [l6( . Consequently, one 
obtains scaling laws for correlation functions in terms of 
an e-expansion. The knowledge of 0(e) terms in the loop 
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expansion in d < d c yields leading order logarithmic cor- 
rections to the mean field answers in d = d c . 

We will now apply the method described above to 
compute both temporal and spatial scaling exponents of 
P t (N, AV) for A + A -> (d c = 2) reaction. 



A. One-loop renormalization of composite operator 



Dimensional analysis shows that renormalization of the 
correlation function (Yli—i't > ( x i>'t)), where 
quires reaction rate renormalization only |il Il2|. How- 
ever, we are interested in single-point correlation func- 
tions of the form (4> N (x,t)), see (|16fl . The operator <p N 
is called a composite operator for N > 2. It is well known 
that insertion of composite operators under the sign of 
averaging leads to new types of divergences in the corre- 
sponding loop expansion [15|. These divergences cannot 
be eliminated by reaction rate renormalization and re- 
quire multiplicative renormalization of the corresponding 
composite operators. As we will see below, these extra 
divergences are responsible for the multi-scaling of prob- 
abilities P t (N,AV). 

For the theory © the first instance of such a 
divergence occurs in the correlator (4> 2 (y,t)) c = 
lim x —>y((f)(x,t)(f>(y,t)) c as e — * 0. We will see that it is 
specifically this singularity that leads to the multi-scaling 
of Pt(N, AV) for N > 2. The tree- level answer is quoted 
below for d < 2[HEI: 



(0(x,t)0(y,t)) c 



8irg t d 



Vt 



- 1 



0(e° 



(22) 

The averaging (• • -) c denotes connected correlation func- 
tions. It is clear from equation (IL'21) that for x ^ y 
there is no divergence in the limit e — > 0. However if 
x = y, there will be a ^ divergence (as e — > 0) despite 
(1221) being a tree-level answer. This divergence cannot be 
regularized by reaction rate renormalization. It requires 
multiplicative renormalization, where the divergent func- 
tion is multiplied by a renormalizing factor which will 
eliminate the singularity. This procedure violates naive 
dimensional arguments according to which ((f) ) should 
scale as (<f>) ~ f or ^ < 2. Multiplicative renormal- 

ization leads to a non-trivial anomalous dimension which 
depends non-linearly on the order of the composite op- 
erator. Multiplicative renormalization is not required for 
the average density ((f)) which explains its lack of anoma- 
lous scaling 

Let Mjv(i) = {(f> N (t)). Our aim is to renormalize Mn 
at (arbitrary reference) time t . Let g — XtJ 2 be the 
bare (dimensionless) reaction rate. The renormalized 
rate g is also defined at reference time t and is given 
by El: 



l + 9o/g* 



(23) 



where g* = ^f- = 2ne + 0(e 2 ) is the non-trivial stable 
fixed point of the RG flow in the space of effective reac- 
tion rates. The constant d is given by 



C d 



(87r)<V 2 1 2 



(24) 



Cd is regular at d = 2 and takes the value j- . Expressing 
g as a power series in g and substituting into Afjv(t ), 
enables us to obtain M/v as a power series in g: 



M N (t Q ) = n c i(g ,t ) \ c 



N 



71=1 



CN, n g 



(25) 



The expression for Mat(^ ) to 0(g 1 N ) is summarized by 
the diagrams in Fig. The first term in this series is 



-X N 



-X N- 




X N-2 



FIG. 3: The diagrams contributing to (4> N ) to 0(g 1 ' N ) 

just ritf and is given by FigOt. Expanding n c i in g using 
relation (|23|) will generate a term proportional to -^r ~ i . 
This term acts as a counter term to the diagram shown 
in Fig|3t>. It eliminates the singularity in the one-loop 
correction to mean density |12| . The first divergence 
which cannot be eliminated by the renormalized reaction 
rate, appears in the coefficient cjv j.. The origin of this 
singularity is from averaging the composite operator <p 2 . 
Let Zn be the renormalizing factor of Mat. 



Zn = 1 + aN - n 9 n 



(26) 



It is chosen in such a way that the renormalized correla- 
tion function Mn.r is non-singular at e = at all orders. 
By definition Hfl]' 



M NjR (g,t ) = Z N (g,e) ■ M(g,t a ,e) 



(27) 



The coefficients ajv jn are chosen to cancel the singularities 
in the coefficients cj\r jn . This is known as the minimal 
subtraction scheme. In particular 



CtJV.l 



Cam 



(28) 
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where [-] s extracts the singular part of an expression. Let 
us look at the calculation in more detail. Using Fig. © 
we may write 



1 


N 







-U.N (/,/)= — — 1 ~ N ^ 7r£ ^ -9 + 0(g 2 ) + finite 

(_29) ' 

The 0(g) term is computed as follows: we need rv^ 2 \ 
there are mC% diagrams of type c and finally by setting 
x = y in (123) which gives [M 2 (t )]s = Hence 



Z N = 1 



N(N - 1) 
Aire 



■9 + 0{g 2 ) 



(30) 



B. RG computation of Multi-Scaling 



The knowledge of renormalization laws 1|23[) and i|30|) 

can be used to compute Mm, ait) for t > t a as follows. 
The Markov property of the evolution operator U for the 
bare theory, namely U(t,t o )U(t o ,0) — U(t,0), tells us 
that the bare function Mm it) is independent of t Q for 
t > t Q . Hence 

d d 
togj-M(t) = t — [Z N X ■ M N , R (t, t )] = (31) 

The function MM, R (t) is a function of (t, t , g), leading to 

the ansatz: MM, R (t) — to Nd ^ 2 ^ (j-, g(t fj , where $ is a 

function with dimensionless arguments. The choice of t Q 
is arbitrary, but we choose it small enough so that MFT 
is still valid. This motivates the choice of the pre-factor 
to Nd ^ 2 . Upon substitution of this ansatz into l(3~T|) we 
obtain the Callan-Symanzik equation: 



Nd 



M N ,R( t > t o,9)=0 (32) 



7jv(s) 

where the f3 and 77V functions of the theory are given by 

(33) 

flln Z„ NY TV — 11 

7Jv(s) 



P{9) = ~to§t = ~ [C d g 2 - eg] 



9g 



The initial condition (at time t ) is given by the loop 
expansion of ((f) ) with the most dominant contribution 
coming from the MFT answer: 



M N ,R(t ,g) = Ki{t ,g) 



(35) 



Equation 132(1 subject to initial condition H35|l is solved 
using the method of characteristics and has the following 
solution for d < 2: 

ft \ Nd / 2 

M N , R (t,g(t,t )) = ijj n N {t ,g(t,t ))- 



g(t,t )-g* 



9-9 



N(N-l) * 



l(t,to) = 



(36) 
(37) 



The running coupling g(t, t ) is the effective RG flow of 
the reaction rate. It can be verified that lim^oo g(t,t a ) = 
g* , which is of order e. Hence for large times we can 
convert the loop expansion to an e-expansion. To obtain 
answers in d = 2 we take the limit e — > in equations 
(OH, ED: 



M N , R {t,g(t,t )) 



g(t,t ) = 



For large times g(t,t ) 



N 



n N {t ,g{t,t )) 



g{t,t ) 



-fin 



(38) 
(39) 



4tt 



Then in the large 



time limit we obtain the following scaling in t (recall 
Mm At) = (4> N (t)) R ) : 



t ~Nd/2 t - N(N - 1 ^ +0(e 2 ) 



d<2 



MM, R (t) \ \ 

(40) 

Combining (|14|l and (|40fl we obtain the following results 
for P t (N, AV) 



Pt(N, AV) 
Pt(l,AV) N 



-+C(e 2 ) 



d< 2 



+ d = 2 

where the AV dependence is restored using dimensional 
arguments. The physical interpretation of (|41() is that in 
the large-time limit particles are anti- correlated (recall 
equation given the same average density the proba- 
bility of finding N reacting particles in AV goes to zero 
faster than the probability of finding N non-reacting par- 
ticles in AV. The origin of anti-correlation can be traced 
back to the recurrence property of random walks in d < 2. 

Let us compare (|41() with exact results for N = 1,2. 
For N = 1 there is no anomaly and our formula is in 
agreement with the well-known result found in |12| for 
d < 2. Let us examine the case N = 2 by studying 
the Langevin SDE : dt<fi — — 2A</> 2 + Ito noise term 
which is equivalent to the field theory (JSJ. Tak- 
ing averages of both sides yields the first Hopf equation: 
d t (<j)) = -2\(<f?). For d < 2 we know (0) - t~ d i 2 . 
Substituting this into the Hopf equation gives ((f> 2 ) ~ 
t -(i+d/2) = t -(d+e/2)_ In d = 2, (</>) ~ l2|i and using 

this in the Hopf equation yields ((f) 2 ) ~ These are 
exact relations and are in agreement with our formula 
gUJ) with 0(e 2 ) = 0. We conjecture that for N = 2, 
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the C(e 2 ) corrections are absent in l|40[) . Finally, note 
from equation 1|23[) lim^oo g = g* . Thus for d < 2, the 
limit A — > oo for finite t yields the same asymptotics as 
the limit t — > oo for finite A. Therefore the large time 
asymptotics of the model at hand belongs to the univer- 
sality class of instantaneous annihilation. 



Comparison of results of £-expansion with exact 
results in d=l 



Equation (|41ll gives us an asymptotically exact result 
in d — 2 (e = 0). In the previous section we have also 
confirmed that order-e expansion of scaling exponent of 
P{2, AV) given by the first equation in 141|) is exact in 
all dimensions d < 2. As it turns out, (|41|l yields an ex- 
act answer for multi-scaling of probabilities in d = 1 for 
all values of N. We are fortunate to have some exact re- 
sults for multi-point correlation functions in the problem 
of diffusion- limited annihilation A + A —> ® for d = 1, 
A = oo [fj . For more details please refer to [(J 0, • Us- 
ing Mathematica and recurrence relations derived in [|| 
we were able to compute exact analytical expressions for 
P t (TV, AV) for N = 1, 2, 3, 4. Based on the results wc find 
that the 0(e 2 ) terms are absent in the e-expansion l)4ip. 
This leads us to conjecture that in d = 1 the one-loop 
answer for the scaling exponents are exact. 

To substantiate our claims, let us review the key re- 
sults from The reaction rate is infinite, hence we do 
not expect to find more than oneparticle at a given site. 
We will use the notation used in The correlation 

function Pn(%i, ■ ■ ■ , %n] t) represents the joint probabil- 
ity density of finding N particles positioned at x\, . . . , xn 
at time t. In particular, p(t) = pi(t) is the average den- 
sity. There is also the convention that Xi < Xj for i < j. 
In the limit of large times: 



P t (N,AV) = dxi...dx N p N {xi 
Jav 



,x N ;t) (42) 



In the large time limit the following answers hold true: 
1 



p(t) 

P2(x 1 ,x 2 ;t) 
P 2 (t) 



(43) 



= 1 - e~ 22 " + V7rz2ie"^ierfc(z 2 i)( 44 ) 

3j t 3j ■} 



(45) 



where D is the diffusion coefficient . Note that as x 2 — > 
x\ (or vice- versa), the correlation function vanishes. This 
is a reflection of anti-correlations between particles. For 
small separations, it can be easily shown that 



p 2 {xi,x 2 ;t) 
P 2 (t) 



(46) 



The above expression for p 2 is valid in the limit of large 
time and fixed separation A = x% — x\, as in this limit 



Z21 — > 0. Due to H42|) the scaling for (|46|l agrees with our 
answer for N = 2 in (|41|) : 



Pt(2,AV) 

p t {i,Avy 



A 

J1J2 



(47) 



With the aid of Mathematica and using similar argu- 
ments as above, it can be shown that P 3 ( Xl f 2 ' X3 ' t '> ^ t~ 3 / 2 

p 3 (t) 

and P 4 ( Zl ^^3,g4,t) ^ f-3 ( which are in agreement with 
our conjecture for the cases N — 3 and N = 4 respec- 
tively. 

The formal reason for anomalous scaling for N = 2 
is vanishing of the two-particle probability distribution 
function at x\ = x 2 . This is the most clear indica- 
tion of anti-correlation between annihilating particles. 
This same phenomenon is responsible for zeros of multi- 
particle distribution functions as well. To explore the 
nature of these zeros starting from the exact recurrence 
relations of we were forced to use Mathematica: exact 
expressions for correlation functions are written as linear 
combination of a variety of terms involving products of 
exp(— z 2 j), erfc(zji), Zji. At small separations these ex- 
pressions simplify due to a number of cancellations which 
are not obvious. Using Mathematica to break through 
the tedious computations we find: 



P2{xi,x 2 ;t) 

P 2 (t) 
P3(x 1 ,x 2 ,x 3 ;t) 

P 3 (t) 
p 4 ,(x 1 ,x 2 ,X3,x 4: ;t) 

P 4 (i) 



= V7TZ21 (48) 

= &Vttz 2 iz 3 iz 32 (49) 
= 2irz 2 iz 31 z 32 Z4iz i2 z 43 (50) 



Note that all distribution functions above vanish as the 
first power of separation between any pair of particles. 
Hence we can make a conjecture about spatio-temporal 
behaviour of distribution functions for arbitrary N: 



Pn(xi, ...,XN]t) 

P N (t) 



n 



(51) 



l<i<j<N 



The above expression is a direct generalization of l|5U|) 
based on permutation symmetry and self-similarity. It 
states that the spatial dependence of the probability den- 
sity is given by the Van-der-Monde determinant of par- 
ticles' coordinates. 

Conjecture (|51|l reproduces the temporal scaling de- 
rived in l|4(J|) as the right hand side contains N(N — l)/2 
factors Zji each contributing i -1 / 2 to the scaling. We will 
present a rigorous proof of this conjecture in a separate 
publication. 



IV. REACTION 3A -> IN CRITICAL 
DIMENSION d = 1 

In the final section of this paper, we extend the pre- 
ceding analysis to the problem 3 A — > 0, where reactions 



occur in triples. Unlike the binary case we cannot ap- 
ply EIM to analyze this reaction in d = 1. Thus we will 
adopt the field-theoretic approach. We give a brief in- 
troduction to the important quantities in the model |12| . 
The action S given by 



S 



d d x / dt 4>{d t - A)0 + 3A00 J + 3A0V + 



\4> 6 (t) A - n <)>5(t) 



(52) 



The Feynman rules are shown in Fig 01 The action must 
be dimensionless which requires the following: 



[t]=L 2 [<f>] 



-d 



[$\=L [A] 



-2d-2 



(53) 



(x,t) (x'.t") 

G (x-x', t-f) 





FIG. 4: Feynman rules for 3A 



-field theory 



The reaction rate is dimensionless at d = 1, which is 
the critical dimension for this problem. Hence in d = 1 we 
expect the system to be characterized by MFT with loga- 
rithmic corrections. As for the binary system, we use the 
classical (MF) versions of the density and response func- 
tions to study the loop expansion. The integral equations 
satisfied by the classical density and classical response 
function are shown in Fig[5| Their solutions are given by 



n c i(t) 

G c i(x 2 ,t 2 ;x 1 ,ti) 



1 



1/2 



6 At 

n c i(t 2 ) 



n c i(h) 



(54) 



G (x2,t 2 ;xi,t 1 )(55) 



Power counting shows that reaction rate renormaliza- 
tion takes care of singularities of Greens' functions which 
do not contain any composite operators. Moreover it 
turns out that all divergences of (4> 2 ) are also eliminated 
by reaction rate renormalization. For higher order com- 
posite operators one needs multiplicative renormaliza- 
tion. 

Let g = Xt e Q be the bare (dimensionless) reaction rate, 
where t Q is our arbitrary reference time. Note that here 
e = 1 — d. To eliminate divergences associated with reac- 
tion rate renormalization we need an expression for g a in 
terms of g. The relation ^ holds, but g* = ^+0{e 2 ). 



a) 



b) 



< 



< 



FIG. 5: a) The classical density; b) The classical response 
function 



The lowest order diagrams contributing to M/v = {<f> ) 
are shown in Fig El The only divergence at the 1-loop 
level associated with composite operators comes from 
the connected three-point function. This determines Zn 
and hence the anomalous scaling. The Callan-Symanzik 



-X 2 



-X N 






-X N-2 



-X N-l 




X N-2 




X N-3 



FIG. 6: The diagrams contributing to ((f> N ) to C^g 1 -^ 2 ) 



equation in d = 1 takes the form 



where 



27T'' 



M NtR {t,t o ,g) = (56) 



(57) 



N(N-l)(N-2) 2 
lN{g) = g + 0(g ) (58) 

Once again we choose the reference time t Q small so that 



9 



the dominant contribution at t Q is given by MFT 

M NtR {t ,g) = n%(t 0) g) 
The solution of l|56[) is given by 



(59) 



t \ N/2 

M N , R (t,g(t,t )) = [jj n%(t 0) g(t,to)) 



g(t,t ) 



JV(JV-l)(JV-2) 



g(t,t ) 



i + m f ± 



(60) 
(61) 



Recall Mjv,ii(t) = {(j) N (t))n. For large times, we obtain 
the following scaling behaviour of composite operators: 



(<p N (t))R 



A 



N 



hit 



N/2 



Out) 



-N(N -1)(N -2) 



1 + 



where 



1 



N 




3N -N(N-l)(N — 2) 



(62) 



(63) 



Using the relation (|14|l between composite operators and 
probabilities 



Pt(N, AV) 
Pt(l,AV) N 



In 



AV 



-JVTjV— 1UJV— 2) 




The result (|64|l reflects the anti-correlation of particles 
as stated in equation J2J- We know from 12} that (0) ~ 

1 /2 

which agrees with l|62|l upon setting N = 1. 
There is no anti-correlation between pairs of particles as 
reactions only occur in triples. Therefore the anomaly 
should vanish for N = 2 which agrees with l|62|) . i|64[l . 
Due to the absence of singularities associated with (<j) 2 ) 
we know the dominant contribution comes from the dis- 
connected n c i diagrams for N = 2. Then (0 2 ) ~ ^ 
which is also in agreement with l|62[) . Finally to check 
(|62|l for N = 3, we can use the first Hopf equation for 

the theory: d t {4>) = -3A(0 3 ). As ((f>) - [ ]s r] 1/ \ the 



Hopf equation implies that (4 
This agrees with (|62|) as well. 



t" 3/2 (lni) 1/2 exactly. 



scaling exponent in the binary system reflects pairwise 
anti-correlation. For the ternary case, the scaling expo- 
nent is cubic which shows anti-correlation within parti- 
cle triples. We obtained our results in a field-theoretic 
setting by identifying probability distributions of multi- 
particle configurations, at scales much smaller than corre- 
lation length, with composite operators in Doi-Zeldovich 
field theory. The origin of the multi-scaling can there- 
fore be traced back to the anomalous dimensions of the 
corresponding composite operators. 



We obtained exact logarithmic corrections to scaling 
for the binary system in d = 2 and for the ternary reac- 
tion in d = 1. We computed scaling exponents for the 
binary system in d < 2 using e-expansion. By analyzing 
the first Hopf equation for the binary system, we proved 
that the one-loop e-expansion gives the exact answer for 
the probability of finding two particles in a fixed volume. 
A similar computation for the ternary system confirms 
the result of the RG computation for the probability of 
finding three particles in the fixed volume. 



RG analysis led us to several conjectures for scaling 
exponents for the A + A — > system in d = 1. First, 
by comparing one-loop answers with the exact results in 
one dimension [6| for TV = 1,2,3,4, we conjecture that 
two and higher loop corrections are absent in d = 1 for an 
arbitrary number of particles TV, see ljH)l. Second, basing 
on the previous conjecture and dimensional analysis we 
propose that the spatial dependence of the multi-particle 
probability density is given by the absolute value of the 
Van-der-Monde determinant of particles' positions, see 
1)51(1. We have recently found a rigorous proof of the 
stated conjectures, which will be published separately. 
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